Entropy production in diffusion-reaction systems: The reactive random Lorentz gas 
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We report the study of a random Lorentz gas with a reaction of isomerization A B between 
the two colors of moving particles elastically bouncing on hard disks. The reaction occurs when 
the moving particles collide on catalytic disks which constitute a fraction of all the disks. Under 
the dilute-gas conditions, the reaction-diffusion process is ruled by two coupled Boltzmann-Lorentz 
equations for the distribution functions of the colors. The macroscopic reaction-diffusion equations 
with cross-diffusion terms induced by the chemical reaction are derived from the kinetic equations. 
We use a if-theorem of the kinetic theory in order to derive a macroscopic entropy depending on 
the gradients of color densities and which has a non-negative entropy production in agreement with 
the second law of thermodynamics. 
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I. INTRODUCTION 

o 

m 

During the last decades, irreversibilityand transport properties have been intensively studied in low-dimensional 
dynamical systems [ElllallllllllEli I IISIIllIllIllIllIllIl- It has been shown that the macroscopic 
£^ • transport coefficients can be related to the characteristic quantities of the microscopic dynamics Moreover, 
results have also been obtained for the entropy balance in dynamical systems sustaining transport processes of diffusion 
0,0,0], electric conduction @, cross effects 0,0, and shear viscosity 0,0,0- For these processes, the 
irreversible entropy production of noncquilibrium thermodynamics could be derived from the underlying microscopic 
dynamics. 

Beside transport, reaction processes have also been investigated from the viewpoint of dynamical systems theory. 
Reaction-diffusion processes play a crucial role in physico-chemical systems far from the thermodynamic equilibrium 
and they have been much studied at the macroscopic level of description on the basis of the phenomenological 
nonequilibrium thermodynamics |l9j |. 

Recently, several microscopic models of reaction-diffusion processes have been introduced and analyzed in order to 
understand the foundations of the phenomen olog ical assumptions. Multibaker and Lorentz gas models of reaction- 
t-H \ diffusion processes were introduced in Refs. [20L l2ll |22|. In this early version, the reaction was an isomerization 
A ^ B with unit probability upon the passage of the particle on a catalytic side where the reaction occurs. In 
Ref. [23, a spatially periodic reactive Lorentz gas was introduced in which the isomerization A ^± B occurs with 
a probability < po < I when the particle collides on catalytic disks. The catalytic disks are few among the disks 
which compose the Lorentz gas. Spatially periodic Lorentz gas and multibaker models with reaction A ^ B have 
been studied in detail and their diffusive and reactive modes were constructed together with their dispersion relations 
HE • I n this wa y ; the macroscopic reaction-diffusion equations could be derived in a systematic way from the 
, underlying dynamics. The analysis revealed the existence of cross-diffusion effects induced by the reaction. Such cross 
diffusion is typically overlooked in the phenomenological approach [27ll2S| . Moreover, the entropy production of the 
standard phenomenological entropy may be negative because of the cross-diffusion effects. Although this problem 
only happens for extreme particle densities it nevertheless sheds some doubts on the phenomenological assumptions. 

The purpose of the present paper is to clarify these issues and understand the problem of entropy production in 
such reaction-diffusion systems. For this purpose, we consider a random Lorentz gas with a fraction of catalytic disks 
where the isomerization A ^ B occurs with a given probability pq. The color A or B carried by the moving particle 
may correspond to the spin of the particle, in which case the catalytic disks model some spin-flipping impurities in the 
system. For dilute Lorentz gases, we can use kinetic theory and linear Boltzmann-Lorentz equations 29, 30]. Thanks 
to such master equations, we can derive a iJ-theorem which allows us to obtain an expression for the entropy. We 
prove that the corresponding entropy production is non-negative with respect to the time evolution induced by the 
macroscopic reaction-diffusion equations. 

The plan of the paper is the following. The kinetic equations of the reactive Lorentz gas are introduced in Sec. [n] 
where we prove the ii-theorem. The macroscopic reaction-diffusion equations are derived in Sec. 11111 The entropy 
balance is obtained and discussed in Sec. IIVI Conclusions are drawn in Sec. [V] 
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II. THE REACTIVE LORENTZ GAS 



In the random Lorentz gas, fixed circular scatterers correspond to heavy particles and moving point particles to the 
light ones. The system is defined by the density of the scatterers n, their radius a, and by the velocity of the light 
particles v = u(cos<p, siny), the system being two dimensional. Because the magnitude of the velocities is fixed, the 
velocity vector can be characterized solely by the angle p. Here, we consider a random distribution of the disks with 
low density n -C a~ 2 . 

The reactive Lorentz gas [24| consists of two different types of light particles A and B, that have a free flight 
between the collisions. Some of the fixed scatterers act as catalysts, i.e., if an A particle collides with such a scatterer 
it becomes a B with the probability po, and vice versa. The catalysts have a density n r . The reaction scheme is 



A + disk 
B + disk 
A + catalyst 



A + disk , 
B + disk , 
B + catalyst , 



with probability po 



(1) 
(2) 
(3) 



The colors can be considered as the two components of a spin one-half carried by the moving particles. The evolution 
of the system can be characterized by the two distribution functions of the two components f A (f, v) and f B if, v) . 
The integral over the velocities of these functions gives the number density at point r of each component i or B. 
Having in view that the time evolution of the distribution function f A (resp. fg) is also influenced by the presence 
of the other component B (resp. A) that may collide with a catalyst, we obtain the system of equations 



d t f A + v-Vf A = 



d t f B + v-Vf B 



av(n — pon r ) 



2 

avpon r 



dp' 



sin ■ 



Lp -if 



[f A (f,p')-f A (r,p)} 



+ 7T 



dtp' 



sm ■ 



<P~<P 



[f B (f,p')-f A (r,p)} 



av(n — pon r ) 



2 

avpon r 



dp' 



sm ■ 



Lp-p 



[fB(r,p')~f B (r,p)} 



dip' 



sm ■ 



<P~<P 



[f A (r,p')~f B (r,p)} 



(4) 



(5) 



One can observe that the total system of the moving particles A and B follows a diffusive process that can be 
characterized by the the sum of the distribution functions 



f = fA + h 



(6) 



The time evolution of / is given by a linear Boltzmann equation, also known as the Boltzmann-Lorentz equation 
IHIH: 



dtf + v-Vf 



avn 



+7T 



dp' 



p-p 



[/(rV)-/(r>)] 



(7) 



This equation rules the so-called diffusion sector. In order to have the time evolution of both distribution functions, 
we need a further equation beside Eq. J7J). 
If we introduce the quantity 



9 = /a - f B , 

and take the difference of Eqs. |Q| and J5J, we obtain the equation 



d t g + v ■ X7g 



av(n — 2pon r ) 



+7T 



dp' 



p-p 



[g(r, <p') - g(r, ¥>)] - 4avp n r g(f, p) 



(8) 



(9) 



which rules the reaction sector. One can notice that the last equation describes a decay in time of the function g that 
can be separated from the rest of the solution 



g(f,p,t) = c- ia ^ t h(f,p) 



The equation for h reads 



d t h + v-Vh = 



av(n — 2pon T 



+ 7T 



dp' 



sm ■ 



p-p 



[h(r, p>') - h(r, ip)] 



(10) 
(11) 



which has the same form as Eq. (JJJ. 
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A. //-theorem 

The Boltzmann entropy of the system of Eqs. Q and JSJ) is given by the sum of the entropies of the two components: 



(12) 



where Boltzmann's constant is taken equal to unity, fee = 1. 

Inserting to the integral on the r.h.s. the time evolution Eqs. (@J and iJSJ, we get for the variation of the entropy 



d t s 



-v-Vf A + C A {f A J B ) 



where the collision integral C A (f A , f B ) has the form 

<p-<p 



In 



eh 



r 



v-Vf B +C B (f A J B ) 



C A {f A , Jb) 
with the notations: 



dip' 



av(n-p Q n r ) avp Q n r 

o UA ~ J A) H ^ — [Jb - J A) 



J A = fA(r,tp,t) , 

!'a = fA(?,<p',t) , 

Ib = f B (r,ip,t) , 

S'b = f B (f,<p',t) ■ 



This relation can be written in form of balance |27 

d t s -- 

where 



~V • J s + (T s . 



is the entropy current and 



cr s = — I dipdip' 



■Is 



ip-lp 



f° f° 
dtp v [ f A In — + f B In — 

JA Jb 



(n- Po n r ) Uf A -f A )l n & + (f' B -f B )lnjk 



+p n r 



(f B -f A )\njk + (f A -f B )lnjA 



> 



(13) 



(14) 



(15) 
(16) 
(17) 
(18) 



(19) 



(20) 



(21) 



is the entropy production. One can see on this general form of the entropy production that it is positive definite under 
the consistency condition that 



< Pon r < n . 



(22) 



III. DERIVATION OF THE MACROSCOPIC REACTION-DIFFUSION EQUATIONS 



A convenient way to find the solution of Eqs. (7J and (|llfl 



dtX + v-VX 



dtp' 



■ p-p 
sin 



[X(r,ip')-X(r,v)} , 



is by writing the distribution functions X — f or h as Fourier series in the velocity angle ip 

+ OC 



X(r,ip,t) = x l(r,t)e 



iltp 



(23) 



(24) 



l = — OO 
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As a consequence of Eqs. (J7J) or (|llfl . the Fourier components Xi fulfill the following coupled differential equations 

d x — id v „ d x + id v „ 8/ 2 c 



1 - 



with 



and 



avn 



for * = 



au(n — 2pon r ) 



for X 



(25) 

(26) 
(27) 



The distribution function X being a real function, the Fourier coefficients with negative index are the complex 
conjugates of their positive index counterpart: X-i = X{ for any If we take solutions in the form 



Xi(r,t) = Ui exp(-7*) exp(ig- r) 



we obtain the eigenvalue equations 

-lUi 



iv — -Ui-i 



ix + iq y 



Ui+i = 



2 2 1-4Z 2 

for / € Z. Expanding in powers of the wavenumbcr q, the decay rates are given by 

8l 2 c 3{Al 2 -l)v 2 2 



Ui , 



7 



4/ 2 - I 



16c 



(28) 



(29) 



(30) 



with / = 0,1,2,... and q = \Jlp '. Figure Q shows the spectrum obtained by solving numerically the eigenvalue 
equations l|29|) for c > 0. We observe that the branch for I — has a convexity which is opposite to the one of the 
other branches for I = 1,2,... Moreover, all the branches terminate on the line 7 = 2c, which fixes the maximum 
value of the wavenumber for each branch. This feature prevents the eigenvalues to become of the opposite sign hence 
avoiding instability in agreement with the stability provided by the _ff-theorem. For the diffusion sector, the constant c 
is always positive. However, the sign of the constant c may change in the reaction sector and we must treat separately 
the cases c > and c < 0. 



' 1 


1 1 1 






1=2 7 




\ --"1=1 




1 


1 



q 



FIG. 1: Eigenvalues -7 of Eq. fgfr versus the wavenumber q = q x for q y = and c = 1. 
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A. The case < pon r < ~ 
In this case, we have that c > in both the diffusion and reaction sectors. 

After the relaxation time, the dynamics is dominated by the first Fourier components Xq, X+i, and X—\: 

X ~ X + X +1 e +lv + X_i e~ ilp . (31) 
In view of the latter relation, we introduce the densities 



P+ = PA + Pb = j f(r, tp) dip = 2irf , (32) 

J — 7T 

P- = PA- Pb = / g(r, V) d<p = 2ng Q . (33) 



Thus the zeroth-order expansion of the distribution function / can be related to the total density. Similarly, we 
introduce the currents 



J+ = Ja + Jb= / vf{r,ip)dip, (34) 

J — IT 
p + Tt 

J- = J A - Jb = V g(r, ip) dip . (35) 



/ * ^\P + + - 2 V-J + ) , (36) 



The first-order Fourier components X + i and A_i are related to the current. As a consequence, we obtain the 
distribution functions / and g in terms of the corresponding densities and currents 

1 / 2 
2n \ v z 

9* JL(p_ + . (37) 

The distribution functions for the species A and B are thus given by 

f A * ^(pa + ^v-Ja) , (38) 

fa * t{ pB + ^- fB ) ■ m 
Equations (|25() for I = and I = ±1 then lead to the coupled equations: 

dt P+ = -V ■ J+ , (40) 

dtJ + = -^nJ + -^V P+ , (41) 

d t p- = —V • J_ — 4avp n r p- , (42) 

*^-?(»-^/-^-. («, 

We notice that the currents relax on a fast time scale so that we can assume that they quickly adjust to their value 
in a quasi-stationary state as 

J+ = -DV P+ , for t » , (44) 

8avn 

J- = -D^, fo r t»^ ¥J , (45) 



and 



with the diffusion coefficient 
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and the reactive diffusion coefficient 

£>« = ^ (47) 

Equation 144(1 is the expression of Fick's law for the particles while Eq. (|45() is its reactive analogue. Substituting 
Eqs. (|44l) and (|45|l into Eqs. (|40|) and 142|) for the densities, we obtain the diffusion equation 

d tP+ = DV 2 p+ , (48) 

as well as the reaction-diffusion equation 

d t p_ = D^V 2 p- - 2k P _ , (49) 

with the reaction rate constant 

k = 2avpan r . (50) 

This result shows that the reaction rate constant is the product of the speed v with the cross section 2a of the disks, 
multiplied by the density n r of the catalytic scatterers and weighted by the probability of reaction p$. The two 
equations (|48|l and l|49[) show the existence of two slow modes in the system corresponding to the decay rate l|30[l with 
I = 0, namely, the diffusive mode of dispersion relation 

diffusive mode: T = Dq 2 + 0(q 4 ) , (51) 

and the reactive mode of dispersion relation 

reactive mode: T = 2n + D {r) q 2 + 0(q 4 ) . (52) 

The equations of motion (|48|l and l|49|) determines the time evolution of the densities pa = (p+ + P-)/2 and 
Pb = {p+ — P-)/2 according to the coupled reaction-diffusion equations: 

d t pA = D AA V 2 p A + D AB V 2 p B - k{pa- Pb) , (53) 

dtps = D BA V 2 p A + D BB y 2 p B + h(pa - Pb) , (54) 
where the transport coefficients can be identified as 

D + Z)W , x 

Daa = D BB = , (55) 

D - L>W , x 

Dab = D BA = . (56) 

The important conclusion is here that there appears a phenomenon of cross diffusion which is induced by the reaction. 
Indeed, the cross-diffusion coefficient takes the value 

D ^=D BA = - Uan{n _ B ^ y (57) 

which vanishes with the reaction probability p a and the density of catalysts n r . In the absence of reaction, the 
cross-diffusion terms vanish with the reaction term and we recover two uncoupled diffusion equations for A and B 
particles. 

We notice that the coupled reaction-diffusion equations l|53|) - (|54|l can be rewritten as 

dtP A = -V • J A - n(p A - p B ) , (58) 

dtPB = -V • J B + n(p A - p B ) , (59) 

in terms of the currents: 

J A = ^ VpA 2 PS ' ( } 

D + D - £>W - , s 

Jb = ^ VpB 2 PA - ( } 

Besides the diffusive and reactive slowest modes, we also find faster modes often referred to as kinetic modes. All 
these modes exist in both the diffusion sector ruled by Eq. (|23f) for X — f and in the reaction sector ruled by Eq. 
(I23|l for X = h. All these modes are characterized by dispersion relations which form a whole spectrum. Figure [21 
depicts the whole spectra of the diffusionand reaction sectors in the case < pon r < ^. 
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(b) 



-r 



FIG. 2: Eigenvalues — T versus the wavenumber q = q x for q y — in the regime < pon r < ^ where c > 0: (a) spectrum of 
the diffusion sector; (b) spectrum of the reaction sector. The parameter values are av = 1, n = 1, and pon r = 0.125. 



B. The case pon r = ~ 

In this case, we remark that Eq. (12;i[) for X = /i has a vanishing coefficient c = in the reaction sector so that the 
equation for X = h is purely advective 



dth + v- Vh = . 
Its solutions are given by h — h(f— vt) and g — e~ 2Kt h(r — vt) obeys 

d t g + v- V.g = -2ng . 

In this case, there is no reactive diffusion coefficient which characterizes the reactive process. 



(62) 



(63) 



C. The case ^ < pon r < n 

As we noticed before, pon r cannot exceed the value n for consistency. In this case, we have that Eq. I|23|l for X = f 
still has a coefficient c > in the diffusion sector but Eq. (|23|l for X — h has a negative coefficient c < in the 
reaction sector. Accordingly, the spectrum shown in Fig. ^ is upside down in the reaction sector and the slowest 
reactive mode is no longer the same as before. 

Here, we must consider the decay rate (|30|) with 1 = 1. The dispersion relation of the reactive mode is now given 

by 



with the new reaction constant 



and the new reactive diffusion coefficient 



T = 2k' + £> (r) V 



c = — (2ra -p n r ) 



9v 



16a(2po"-r — n) 



(64) 



(65) 



(66) 



A transition should therefore appear at high concentrations of catalysts. 

Figure |21 depicts the spectra in the diffusion and reaction sectors in the case 5 < pon r < n. 

In the following, we shall focus on the case < pon r < \ where the catalysts are dilute enough. 
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IV. THE ENTROPY BALANCE 



In this section, we derive the equation for the balance of entropy from the ii-theorem in two approximations for the 
entropy density. In kinetic theory, we have a well-defined expression for the entropy which guarantees that the entropy 
production is always non-negative. However, at the level of the macroscopic description given by the reaction-diffusion 
equations (|53[) - i|54|l . the entropy is given as an approximation of the expression (|12[l of kinetic theory and we must 
verify the domain of validity where the corresponding entropy production is non-negative. 

The first approximation we consider is based on the standard phenomenological entropy defined in irreversible 
thermodynamics of reaction processes [27], . We show that the corresponding entropy production is non- negative 
in a broad range of values for the densities pA and ps including the thermodynamic equilibrium state but there is a 
small domain where the entropy production corresponding to this approximation fails to remain non-negative. 

Therefore, we consider a second approximation which includes extra terms involving the gradient of the densities. 
We show that the corresponding entropy production always remains non-negative. 



A. Entropy without gradients 



Supposing that the system is sufficiently dilute, the phenomenological irreversible thermodynamics supposes that 
the entropy density has the following expression: 

p° p° 

s = p A In 1- p B In — . (67) 

PA PB 

This entropy density is obtained from the entropy (|12[1 of kinetic theory by using the expansions (|38fl and (|39fl of the 
distribution functions and by keeping the terms in the densities themselves and discarding terms in the gradients. The 
reference density is thus given by p° = 2irf°, which amounts to suppose the equality of the masses of the particles A 
and B. 

The variation of the entropy density s in time is given by 

d t s — V • J s o s , (68) 

as calculated by using the coupled reaction-diffusion equations (|58|I - H59() . The entropy current density is obtained in 
terms of the currents (|BTi|) and l(6"T|) of particles A and B as 
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while the entropy production as 



a s = K(p A - p B ) In J a- Jb 

PB PA PB 

Inverting Eqs. (|60|) and l|tjl|) . we obtain the gradients in terms of the currents as 



J A 2k 
~D + 3^ 
J B 2k 
U ~ 3^ 



VPs = -—-^(Ja-Jb) 



Substituting in the entropy production (|70|) . we get 



k(pa - Pb) In h aJ A + 2j3J A • Js + . 



(70) 

(71) 
(72) 

(73) 

(74) 
(75) 
(76) 

The first term in the entropy production l|73|) is always non-negative because {x — y)hi(x/y) > for all positive 
values of the real numbers x and y. On the other hand, the last terms constitute a quadratic form which is non- 
negative under the conditions that /3 2 — aj < with a > and 7 > 0. For the diffusion coefficient l|4t)[l . we check 
that 



with the coefficients 







2k \ 


PA y 




3t^y 


K 


(± 




3^2 


VPA 


f PB/ 






2k \ 






3^7 



7 



4a 

3wpA 
4a 

3v/5b 



(4n.-p ^r) > , 
(4n-p n r ) > , 



because of the consistency condition (|22|l . Next, the condition /3 2 — < is given by 

2" 



e 2 



1 - 2 



in terms of the density ratio 



and the parameter 



1 



PB_ 

PA 



£, + K 



P0«r 



4k£> 



(77) 
(78) 

(79) 
(80) 
(81) 



2n 3v 2 

The roots £± of Eq. (|79|l satisfy £+ = jT-- If £+ > the entropy production is thus non-negative in the domain 

C+PA > PB > Z-PA ■ (82) 

It turns out that the sum of the roots 



t+ + £- = 2 



21 --1] -1 



(83) 
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is always positive in the interval < y < i where the condition of consistency i|22|) is verified. (The sum £ + + £_ is 
negative only in the interval 4 — 2y/2 = 1.1715... < y < 4 + 2y2 = 6.8284... outside the domain of consistency.) Since 
£ + and have the same sign, they both are positive in the consistency domain < y < |, which means that there 
exists a domain of the quadrant (pa > 0, ps > 0) where the entropy production (|73|) can be negative. This domain 
is composed of £_ pa > Pb > and £_ pb > Pa > 0. Another way of saying it is that the domain lj%2l where the 
entropy production is non-negativity is smaller than the quadrant {pa > 0, pb > 0) of all the possible densities /cu 
and pb- The domain of non-negativity extends to 

1441.9... pa > Pb> 0.00069348... pa , for y = 0.1 , (84) 

to 

321.99... p A > Pb > 0.0031056... p A , for y = 0.2 , (85) 

and to 

33.970... pa > p B > 0.029437... p A , for y = 0.5 . (86) 

It is only in the limit y = without chemical reaction that the domain of non-negativity coincides with the quadrant 
(pA > 0, pb > 0). In the presence of the chemical reaction, the situation is unsatisfactory because the scheme is not 
consistent with the second law of thermodynamics even if the domain of negative entropy production is small and 
only concerns color densities which are very far from the thermodynamic equilibrium. 



B. Entropy with gradients 

To cure the problem reported here above, we choose to expand the entropy production (|12fl of kinetic theory to 
include the terms with gradients of the densities. For this purpose, we substitute the expansions l|38|) and i|39[l of the 
distribution functions in the entropy (|12|l and we truncate up to the terms which are quadratic in the currents to get 

S = PAm^+p s m^-lpU^ , (87) 

PA PB V- ypA PB J 

where the particle currents are given in terms of the density gradients according to Eqs. (|60|l and (|61|l . Accordingly, 
the entropy l|87|l is quadratic in the density gradients. 

Here again, we calculate the time variation of the entropy density s by using the coupled reaction-diffusion equations 
(IHEJl-lEnjl- The balance equation for entropy is here also given by 

d t s = -W-J s + a s , (88) 

with an entropy current similar to Eq. (|69|l 

P° - P° 

J s = J A ^— + J B ln-^—+ 0(3) , (89) 
ep A ep B 

up to the terms 0(3) of third order in the gradients. However, the entropy production now takes the more complicated 
form 

v a = K(p A -p B )lnB±-J A .Y£A-j B .Ym + ^ { j A -j B) .( -^(pa-pb) (4 - 41 +0(4) , (90) 

pb pa pb v z y PA pb J v 1 y PA p B J 

where 0(4) denotes terms of fourth order in the gradients. Replacing the gradients by the currents with Eqs. (|71|) - (|72|l . 
we obtain 

a a = .(pa -pB)\n P -^ + ^( J l + J l\+ ^- 2 (J A - Jb) ^(p A - Pb) ( 4 " 4) + 0(4) , (91) 

p B d y PA p B J 3v 2 yp A p B J v 2 y p\ p B j 

or equivalently 

o s = k( Pa - pb) In P — + a ,P A + 2f3J A ■ J B + <yJ% + 0(4) , (92) 

PB 
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with the coefficients 



1/1 K 

2k ( 1 



3v 2 \pA 
1/1 K 



k Pb \ 
1 + ^— ) > 


(93) 


v z p A J 






(94) 


PbJ 




(95) 


1 v 2 p B ) 





p B \D 3v : 

The coefficients a and 7 are always positive, while the condition 1 — < of non-negativity of the quadratic form 
is here given by 

£ 2 + 2^±|^ + l>0, (96) 
2/(12 -y) 

with the ratio £ defined by Eq. I|8(J|I and the parameter < y < | by Eq. (|8ip. The opposite inequality is obtained 
compared to Eq. i|79|) because the dependence on the ratio £ is here more complicated but still simple enough to lead 
to the quadratic equation (|96|l . In the physical domain < y < | , the roots £ + and = l/£+ of Eq. H96|) are real 
and negative because 

£+ + l _ = _ 2 (£±|>l±l<„. (97) 

2/(12 - y) 

As a consequence, the quadratic part of the entropy production l|92|) is non-negative in the quadrant of all the 
physically allowed densities where pa > and ps > 0. If the gradients of densities are small enough so that the 
corrections of fourth order in Eq. I|92|) are negligible, the whole entropy production (|92fl is also non-negative. 

We have thus proved that the inclusion of the gradient terms in the entropy avoids the aforementioned problem and 
guarantees that the entropy production remains non-negative for all the values of the color densities if the gradients 
are sufficiently small. 



C. Interpretation of the gradient terms in the entropy 

The entropy density given by Eq. (|87() can be expressed as follows in terms of the gradients of the particle densities 
by using Eqs. 

s = p A hi — + pb In — — ^ NpaY - KabVpa ■ V P b - ^ Npb) 2 , (98) 
p A Pb 2 \ / 2 V / 

with the coefficients 

(D + D^) 2 (D-D^) 2 , x 

K AA = '- + ^r- 2 L , (99) 

(D + D^) 2 (D-D^) 2 

2v 2 \p A PbJ 

which are independent of the velocity v. 

The gradient terms are of the same kind as those appearing in the Ginzburg-Landau free energy. Here, they appear 
in the entropy with the opposite sign in agreement with the required thermodynamic stability of the equilibrium state 
|3lL l32j|. Indeed, the entropy must be maximal in a stable equilibrium state which is the case since the quadratic 
form in Eq. (|87|l or (|98|l is negative. Accordingly, the entropy reaches a maximum at the equilibrium state where the 
gradients and the currents vanish. 

The gradient terms are responsible for statistical correlations between the particles. Indeed, the entropy density 
H98JI can be used to define the entropy functional 



S[pa,Pb] = J sdr 



(102) 
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and the probability distribution for statistical average given by the functional integrals 



/ V P aD P b O exp ( j-S[p A , Pb) 
iO) oq = ^ , (103) 



/ Vp A Vp B e^p(-^S\fiA,PB. 



for an observable O. This allows us to calculate the correlation functions of the particle densities at the thermodynamic 
equilibrium. We consider the correlation functions of the densities <|32[) - (I33[) : 



C++(r) = (p + (r) P+ (Q)) cq - (p + )i q , (104) 

C__(0 = (p-(^P-(0)) cq -( P -) 2 cq , (105) 
C+_(r) = ( P+ (r)p_(0)) cq -(p + > oq (p_) eq . (106) 

Partial differential equations can be obtained for these correlation functions by the variational principle SS — based 
on the entropy functional (|98() evaluated around the thermodynamic equilibrium 



(p + -p eq y + pi K ++ , cq ^ K , _, , 2 

2p eq 2 



V P+ -^T^ Vp- , (107) 



where the equilibrium values of the coefficients are given by (p+) C q = Pcq = \p , (P-)eq = 0, s eq = p C q, and 

1 2_D 2 
i^++,c q = 7 {Kaa + K BB + 2K AB ) cq = , (108) 

4 V p cq 

1 2D^ 2 
K— ,eq = 7 (#AA + K B b - 2K AB ) cq = — , (109) 

^+-, cq - J {K A a - K BB ) cq = , (110) 

so that there is no cross term in the gradients of p± in Eq. (|l()7(l . The fluctuations described by Eq. (|103|l are 
Gaussian around the equilibrium since the entropy density H107I) is quadratic. As a consequence, the correlation 
functions obey the equations 



C++(r)=0, (111) 



C— (*0=O, (H2) 



for r ^ and with the correlation lengths 



3\/2 

^+ = yfp cq K ++ ^ cq = — , (113) 

loan 

3%/2 



and C-| (r) = 0. These correlation lengths are of the order of the mean free path of the particles between the 

scatterers of radius a and density n. The correlation functions are given by Bessel functions of zeroth order and they 
behave at long distance as 

C±±(f) ~ ±= cxp (~\ , (115) 



V. CONCLUSIONS 



In this paper, we have studied a reactive random Lorentz gas in which a point particle carrying a color A or B 
(or a spin one-half) bounces among randomly distributed disk scatterers. Some of these disks are catalytic in the 
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sense that the reaction A ^ B occurs with a given probability po upon collision on these catalytic disks. In the case 
of a particle with a spin, the catalytic disks correspond to impurities flipping the spin. Under dilute-gas conditions, 
the time evolution of the distribution functions of finding the particle with a given color or spin at some position r 
with some velocity v are ruled by two coupled Boltzmann-Lorentz kinetic equations which satisfy a if-theorem. The 
ii-quantity defines the entropy at the kinetic level of description. 

The time evolution separates into a diffusion sector for the total distribution function and a reaction sector for the 
difference of the distribution functions of the colors. The spectrum of eigenmodes of both sectors can be constructed in 
detail. The diffusive (resp. reactive) modes are the slowest modes among all the modes of the diffusion (resp. reaction) 
sector, which provides us with the diffusion coefficient of the diffusive mode, as well as the reaction rate constant and 
the reactive diffusion coefficient of the reactive mode. This analysis and the knowledge of these coefficients allow us to 
obtain the macroscopic reaction-diffusion equations. These equations present cross-diffusion terms which are induced 
by the reaction if the reaction probability vn is not vanishing, as in the previous studies of the reactive periodic 
Lorentz gas and multibaker models pel l2ll 124 . 125ll26| . In the reactive random Lorentz gas, our analysis shows that 
a transition happens in the reaction sector between a regime at low concentrations of catalytic disks and reaction 
probabilities and another one at high concentrations and probabilities. 

Using the derivation of the macroscopic reaction-diffusion equations, we have studied the problem of the entropy 
production on the basis of the entropy defined in kinetic theory in terms of the distribution functions and the 
associated ii-theorem. The entropy of kinetic theory can be expanded in powers of the gradients of the densities 
of both colors A and B. At the lowest order of this expansion, the entropy density is simply a function of the 
color densities themselves and coincides with the expression of the phcnomenological nonequilibrium thermodynamics 
[13, The balance equation of this entropy without gradient can be derived from the macroscopic reaction-diffusion 
equations. Because of the cross-diffusion effects induced by the reaction, the resulting entropy production may become 
negative for extreme values of the ratio between the color densities. This is certainly a problem of principle for the 
phcnomenological approach. In order to solve this problem, we have considered the entropy density at the next order 
of the expansion in the gradients of color densities. At this next order, the entropy density contains terms which 
are quadratic in the gradients beside the contribution of the phenomenological entropy. The balance equation of this 
entropy with gradients turns out to have an entropy production which is non-negative for all the color densities and 
for small enough gradients of color densities, in consistency with the second law of thermodynamics. 

The entropy density with gradients is interpreted as the entropic version of the Ginzburg-Landau free energy. The 
addition of gradient terms are shown to be responsible for statistical correlations in the densities of the colors A and 
B over spatial scales of the order of the mean free path of the particle. The consideration of these spatial correlations 
appears to be necessary to get a non-negative entropy production in the presence of chemically induced cross diffusion. 
The inclusion of these gradients in the entropy is justified by the kinetic theory and by the consistency so obtained 
with the second law of thermodynamics. Such an entropy functional with gradients is also justifed by analogy with the 
Ginzburg-Landau free energy |3lll32j . The inclusion of g radients in the entropy requires a departure from the classical 
Onsager-Prigogine nonequilibrium thermodynamics |27l |2SJ . This classical nonequilibrium thermodynamics neglects 
the possible interplay between the diffusion and the reaction, such as the cross diffusion induced by the chemical 
reaction. This effect is expected in systems with a high reaction probability and high concentrations of reactants with 
respect to the inerts species which do not participate to the reaction. We think that the present workk clearly shows 
that this chemically induced cross diffusion is compatible with kinetic theory and the second law of thermodynamics 
and should be an experimentally observable effect. 
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